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A variational treatment for a two-electron quantum dot (the artificial helium atom) is proposed 
which leads to exact answer for the ground state energy. Depending on the magnetic field value the 
singlet-triplet and triplet-triplet transitions of the ground state take place, which are captured in our 
solution. Using the same variational technique we find corrections to wave function and energy and 
the transition field srengths in a realistic dot where electron wave function has a finite extent in the 
direction perpendicular to the x-y plane in which usually 2-D dot electrons are confined. Using the 
variational wave function we show that photoemission cross-section as a function of magnetic field has 
sharp discontinuities, which can be used for experimental verification of the singlet-triplet transitions. 

Quantum dots [1] are little two-dimensional islands of electrons, which are laterally confined by 
an artificial potential. They can be thought of as artificial atoms with the field of nucleus replaced 
by an imposed external potential. The artificial hydogen atom is a single electron in a two dimensional 
circular geometry confined by a harmonic potential. The problem becomes interesting in the presence 
of magnetic field in the perpendicular direction and wave funcions for this case were worked out by 
Fock [2] shortly after the Schroendiger equation was established. The artificial 'helium atom' problem 
was, however taken up more than forty years later. In an extensive numerical work Maksym and 
Chakraborty [3] and Wagner et al [4] found extremely interesting effect of the competition between 
magnetic field and Coulomb repulsion between electrons in two-electron quantum dots. In particular 
these authors found that the ground state can change character as the magnetic field changes, leading 
^Author to whom all correspondence should be addressed 
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to singlet-triplet transitions. Six years later Dineykhan and Nazmitdinov [5] using the formidable 
tools of constructing equivalent hamiltonians in oscillator representations solved the problem for a 
two-electron quantum dot exactly. In this communication we note the fact that in the real helium 
atom problem, excellent results for the ground state are obtained by using a variational principle 
and exploit that to set up a variational calculation for this artificial helium atom. The ground state 
energy turns out to be the same as obtained from the exact treatment. Employing the same variational 
techniques we find the required correction for energy, wave function and the magnetic field strengths at 
which transitions are taking place for a more realistic quantum dot in which elctron wave function has 
finite extent in the third direction. This is done by imposing a more stronger harmonic confinement 
in the z-direction. We also show that the photoemission cross-section calculated using the variational 
wave function has sharp discontinuities at the magnetic field values where the transitions occur which 
has been proposed to probe it experimentally [6]. 

For two electrons confined by a parabolic potential in two-dimensional plane,the hamiltonian is 
(gauge = \By, Ay = -\Bx) 

and if the confinement in the z-direction is considered the hamiltonian is modified by the term 

fi2 rP 1 
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and the Coulomb term becomes | ^-|_'^^^ | . In the H above pi and p2 are the two dimensional position 
vectors of the two electrons (n and r2 are three dimensional position vectors) , oJc = is the cyclotron 
frequency, luq is the frequency of the confining potential in x-y plane {coz is the confining frequency 
in the z-direction and lUz S> coq), m* is the effective mass of the electron in the semiconductor, e is 



,2 , 



the dielectric constant and g is effective Lande factor for the semiconductor and = (wq + ^). We 
transform to the center of mass coordinate pc = \{pi + P2) and the relative coordinate p^ei = {pi — P2) 
and write the hamiltonian in equation (1) as = i?c + H^-eU where Hj-^i = H^^i + Hint and 

He = -^^l + Y^-^nv,J + m*n'pl (3) 

H'rel = -^^lel + ^{-if^^,rJ + \m*^'plel (4) 

1 

Hint = (5) 

47reeo Prei 

The wave function '^{pcPrel) will clearly separate as '(pi{Pc)ip2{Prel) with the energy eigenvalue E 
splitting as E = Ec + Ej-ei where 

H.MPc) = EcMp^c) (6) 

Hreli^2{Prel) = Erel1p2{Prel) (7) 
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He and H^^i axe hamiltonians of single electron quantum dots with the masses of the electrons given 
by 2m* and m*/2 respectively. Consequently the exact answers are known for these parts and we have 



Ec = i2N+ I M I +i)hn 
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^0,; = (2n+ I m I - J ^ (8) 



with the wave functions given by 

i^lNMipc) ■ 



(^)l^lLH(^)xexp[-^-zM</..] (9) 



i>2nm{Prel) - ^ 2H+%a|,r(n+ | m | +1) ^" ^2^^^ 

exp[-^f- - %m<^reA (10) 

Here radial quantum numbers N, n are positive integers and angular momentum quantum numbers 
M, m can take all possible integral values. Length scales are set by an and an and = ^^^^ = 
In addition to these if we consider the Zeeman term then there would be an energy contribution 

2 

Espin = 9PbB J2 Sj = qubBSI^^^i (11) 

i=i 

where is the z-component of the spin operator of the j-th electron and g is the effective Lande factor 
for the semiconductor. When permutation of electrons take place prei —Prei and antisymmetry 
requirement implies for odd m values triplet and for even m values singlet states and Espin = giJ,BE{l — 

(-ID- 

However there is no exact answer for Hint and consequently approximation techniques have to 
be resorted to. Our observation is that the variational principle which one employs to calculate the 
ground state of helium atom should be effective here as well. This is what prompts our trial wave 
function. The center of mass motion docs not enter the picture, the ground state of that being fixed 
by = M = 0. In the absence of the Coulomb interaction, the ground state of the relative motion 
would be given by n = m = 0. It is important to note that in the absence of confinement (i.e.ujQ = 0), 
the energy levels arc independent of m, but in the presence of confining potentials, for a given value 
of n, there is an interesting dependence on the azimuthal quantum number m. Magnetic field tries to 
compress the wave function i.e., to decrease the separation of the electrons where as to minimize the 
Coulomb repulsion electrons would tend to increase their separation. Optimization of these competing 
effects i.e., minimization of the total energy takes place at different m values depending on the strength 
of the magnetic field. Now to maintain the anti-symmetry of the total wave function the spin state 
changes between triplet to singlet corresponding to change in m from odd to even value and vice 
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versa as evident from the Egpin expression. At strong magnetic field due to spin polarization, spin 
states would be triplets; thus ground state ceases to show further singlet- triplet transition and only 
triplet-triplet transitions take place. Consequently, in writing down the ground state wave function, 
we respect this fact by anticipating an n = form for the radial wave function but allowing for an 
azimuthal dependence as exp(— ■im^rei)- -^^ the dominating part of the potential terms in H^^i + Hint 
is the parabolic one, the effect of the Coulomb part can be accounted by introducing the length scale 
as the variational parameter in analogy with the helium atom problem. This inspires the variational 
trial wave function (normalized) 



= V 2H.1.;54(|„,K1) '7''""'°"'''-^ - '''' 
With this trial function, there will be four contributions to Erei, which we write as 



^V.,, I i^om >= ^ 



E2{P,m) = < tporn I y(-«^V0^^J I V'Om > ^ 

E3{p,m) = <Vom| J^|V'Om>= ^(|m|+l) 

E^iP,m) = <^om\j |V'Om>=- 7^^, (13) 

Aireoeprei 47reoe/3v2r(| m | +1) 



This leads to 



^ N /, , xr m*ri2^2 io'^mh f-r(| ;» I +i) 

EM = (i 1 + ^1 - — + 4..„./?^2r(| n. I +1) 

Minimization of this energy with respect to /? leads to 

^4 an 1 r(|m|+l) ^ 

a* V2T{\m\ +2) ^ ^ 



where x = and a* = ^^p^- We solved equation (15) numerically for GaAs bulk conduction band 
electron parameters {g = —0.44, ^ = 0.067, e 13) and the results are shown in fig.l and fig.2. We 

E -\-E 

have taken the typical confinement energy to be 4meV and varied B from to 12T. In fig.l 
vs. B is plotted. The transitions m = 0^m = l, m = l^m = 2 occur respectively at magnetic 
field strengths 1.3T, 6. IT (approximately). The mean square separation of two elctrons is given by 
< V'Om I Prei I V'Om >= 2(| m \ +1)/?^ and root mean square separation d = V <^om\p^^i\i>om> g 
plotted for different m values in fig.2 showing the discontinuities brought about in the inter-electron 
separation by the transitions. 

However for a more realistic quantum dot the finite extent of the electronic wave function in 
the third direction has to be taken into account and we achieve this by imposing another harmonic 
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confinement in the z-direction and employ our variational tools to get the characteristic equation. 
With the Hz part we have the wavefunction modified by 

M^) = I ^ eM-^)Hn{\) (16) 

^J2'^T{nz + l)7rV2A A 

and the energy is modified by the term Ez = {nz + \)fii^z- Separating the hamiltonian again in center of 
mass and relative coordinates and following previous arguments it can be clearly seen that the ground 
state would have quantum numbers AT^ = and n^^' = 0. We introduce the variational parameters 
Pi and (32 respectively replacing an and X^ei and after energy minimization get coupled characteristic 
equations for x = and y = where A^g^ = 2A^ = ^^^^ . The characteristic equations are 

" - ^^^^^A;;;) ril+\m\)'^'^r'^ \m\,-+\m\,l- 2(-) (— ) )] - 1 (17) 

4 _ A^ r(|m | 1+ I ^ I, ^+ I m 1, 1 - 2(^)2(^)2) + 2^ f X 

^ a* ^r(|+|m|) '4' ' ''2 ' ^y' ^Xrel y a*XreiT{l+\m\) 

2i^i(|, 2+ I m 1, 1+ I m 1, 1 - 2(^)2(|^)2) = 1 (18) 
2 2 y Arel 

We have solved these coupled equations simultaneously with GaAs parameters and found that singlet- 
triplet transitions are taking place at higher magnetic field strengths and these results are shown 
in Fig.3. These results can be understood from the physical ground that with the introduction of 
harmonic confinement in the z-direction x values decrease significantly compared to the exact 2-D 
situation (Fig. 4) leading to increase in Coulomb energy and to minimize the total energy the required 
transitions need higher energy contribution from the magnetic field. Therefore, the transitions occur 
at higher fields strengths. All these results are in good agreement with the exact solution of Dineykhan 
and Nazmitdinov [5]. 

Usually the energy levels in the semiconductor quantum dots arc probed by the FIR spectroscopy. 
Due to the long wavelength of the electromagnetic field (compared to the dot size), there is no appre- 
ciable change in the electric field across the dot and hence, the elctric field couples only to He through 
the contribution eE ■ (pi + P2) cxp(iujt) = 2eE ■ pcexp{iujt), where £^ is a constant electric field (dipole 
approximation). As Hrei is not coupled to the field, FIR spectroscopy does not probe the effects of 
inter-electron repulsion. However if wc consider the photoemission from the dot to the vacuum by 
irradiating the dot with high energy photons, then Hrei is coupled to the external fields (as the wave 
length is comparable to the dot size, dipole approximation is not valid) and the singlet-triplet and 
triplet-triplet transitions arc manifested as sharp discontinuities in the photoemission cross-section. 
The calculation of the photoemission cross-section involves the evaluation of a matrix element cor- 
responding to the electromagnetic field operator causing the emission from the initial state to final 
state and thus depends on the details of the initial and final state wave functions. This dependence 
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on the initial and final state introduces the m dependence of the cross-section and in this case we do 
not replace ex.p{ik ■ r) by unity as done in the dipole approximation. The operator is treated exactly 
and the final form of the differntial cross-section calculated using our analytically found trial wave 
function becomes 

^ 2S+3m • 2. , Jx2 + 2)2H(i^Q^)2H ^2a|(3x2 + 2) 
— - = 1 — 7T 1 — ( "m +1) — \-— sm 6'cos 6 ^, , o , , „ exp — 2 ^ — jrr — -\ 

(19) 

Here, q is the wave vector of the emitted electron, k is the wave vector of the incident photon, 
hK = h(k — q) is the momentum transferred and 6 and 4> are respectively the angles q makes with 
k and kn plane where h is the unit polarization vector of the incident photon, = ^fjij(^)^ is a 
constant extracted to express the differntial cross-section expression in a dimensionless form. So, 

K'^ = k^ + q^ -2kqcose (20) 

holds and coscj) = sm9nCos{(f)fi — (j)^) and cos 6 = sin cos (0^ — (j)^). Therefore, it is very clear 
from the expression for the differential cross-section that it depends significantly on the direction of 
incidence and polarization. This for some simple cases can be illustrated easily. If k is parallel to 
z-axis then cos 6 = and cos (p = cos (f)^ or cos 4> = sin (pq as h is parallel to x or y-axis. So, the 



angular distribution is proportional to cos^ 4>ff or sin^ cp^ and if it is the case of circular polarization 



then the angular distribution is proportional to (cos (j)n cos 4>^+ sin (j)f^ sin (p^) and only if n = ^ 
then it becomes isotropic. But, when the k lies in the x-y plane then with all the cases of circular 
polarization we shall have angular dependence. It also becomes apparent from the expression that 
emission count is larger in the direction of polarization compared to other cases and if the photon is 
linearly polarized in the z-direction then there is no emission. So, depending on the 'm' values of the 
ground state as a function of magnetic field the cross-section would have different angular distribution 
as well as discontinuities characterizing transitions of the ground state. Based on these one can probe 
now these transitions experimentally. In the context of experiment the frequency of the photon has to 
be carefully chosen as well as those directions mentioned above and details of these shall be provided 
elsewhere [6] where we also show the modifications in the context of more realistic quantum dots due 
to finite thickness. For illustration we plot the above expression for the transition m = — > m = 1 
and m = 1 ^ m = 2 in fig. 5. 

It should be noted that, using the proper variational wave function for the ground state we 
have reproduced the exact solution of Dineykhan and Nazmitdinov [5]. The essential features of 
the interacting two electron ground state of the semiconductor quantum dots, i.e., the singlet-triplet 
transitions with the change in the magnetic field strength and triplet-triplet transitions at very strong 
magnetic field are explained using our result. In the triplet state the angular momentum of the ground 
state is m = 2p + 1, where p is a positive integer and this is weak version of the magic number concept. 
As, in our choice of the trial wave function the symmetry under permutation is properly taken into 
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account and the orthonormalization is built in, proceeding in the same way it is quite trivial to find 
the higher excited states with nonzero radial excitations i.e., n 7^ . The competing effects of the 
magnetic field and Coulomb repulsion are also demonstrated in terms of the variational length scale, 
which captures the discontinuous changes in the inter-electron separation at transition points. Also, by 
same kind of variational calculation the wave function and energy found for a dot with finite thickness 
are in good agreement with the physically anticipated results. Prom the calculated photoemission 
cross section using the above mentioned trial wave function we have also analytically demonstrated 
the sharp discontinuities and the difference in the angular distribution at transition points which can 
be used to probe these transitions for the parabolic dot in addition to the found discontinuities from 
the magnetization measurement. 
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Figure I: Z = vs. B(T) is plotted for different 'm' values and m = 0^m = l, m = l^ 

m = 2 transitions are taking place at i? = l.ST and B = Q.IT respectively and also other energy level 
crossings are present. 
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Figure 3: E = ^''"''y^^^^^" vs. B(T) plotted for different 'm' values when finite thickness of the dot is 
taken into account. For this purpose ^ = 9 has been taken. 
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Magnetic field B 
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Figure 4: x vs. B(T) is plotted for m = for both the 2-D and reahstic 3-D dot and it is clearly found 
that throughout the range of magnetic field x values have decreased for the later one. 
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Magnetic field B 



Figure 5: vs. B(T) is plotted for different 'm' values of the ground state as the B is varied and 
showing discontinuities as characteristic of the transitions. Prom the plot it is found that percentage 
change for m = l— >m = 2is~ 4.5 times smaller compared tom = 0^m = l transition and also it 
is evident from the plot that at a fixed frequency behaviors of different 'm' cross-sections are going to 
change and for this reason both lo and B have to be varied to observe all the transitions properly. 



12 



